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Energy transfer in binary collisions of two gyrating charged particles in a 
magnetic field 

H. B. Nersisyar0 and G. Zwicknage0 

Institut fur Theoretische Physik II, Erlangen-Niirnberg Universitat, Staudtstrasse 7, D-91058 Eriangen, 
Germany 
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Binary collisions of the gyrating charged particles in an external magnetic field are considered within a classical 
second-order perturbation theory, i.e., up to contributions which are quadratic in the binary interaction, 
starting from the unperturbed helical motion of the particles. The calculations are done with the help of a 
binary collisions treatment which is valid for any strength of the magnetic field and involves all harmonics 
of the particles cyclotron motion. The energy transfer is explicitly calculated for a regularized and screened 
potential which is both of finite range and nonsingular at the origin. The validity of the perturbation 
treatment is evaluated by comparing with classical trajectory Monte Carlo (CTMC) calculations which also 
allow to investigate the strong collisions with large energy and velocity transfer at low velocities. For large 
initial velocities on the other hand, only small velocity transfers occur. There the nonperturbative numerical 
CTMC results agree excellently with the predictions of the perturbative treatment. 



I. INTRODUCTION 

In the presence of an external magnetic field B the 
problem of two charged particles cannot be solved in a 
closed form as the relative motion and the motion of the 
center of mass are coupled to each other. Therefore no 
theory exists for a solution of this problem that is uni- 
formly valid for any strength of the magnetic field and the 
Coulomb force between the particles. The energy loss of 
ion beams and the related processes in a magnetized plas- 
mas which are important in many areas of physics such as 
transport, heating, magnetic confinement of thermonu- 
clear plasmas and astrophysics are examples of physi- 
cal situations where this problem arises. This topic was 
studied starting with the classic papers in Ref. ll| using 
kinetic equation approach, where the binary collisions of 
the particles are masked due to the velocity average in 
the collision operator. Recent applications are the cool- 
ing of heavy ion beams by electrons** and the energy 
transfer for heavy-ion inertial confinement fusion.— 

Calculations have been performed for binary collisions 
(BC) between magnetized electrons^ and for collisions 
between magnetized electrons and ions i^r— In the pres- 
ence of a magnetic field only the total energy E of the 
interacting particles is conserved but not the relative and 
center of mass energies separately. In addition, the rota- 
tional symmetry of the system is broken, and as a con- 
sequence only the component of the angular momentum 
parallel to the magnetic field is a constant of motion. The 
apparently simple problem of charged particle interaction 
in a magnetic field is in fact a problem of considerable 
complexity and the additional degree of freedom of the 
cyclotron orbital motion produces a chaotic system with 
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b 'Electronic mail: guenter.zwicknagel@physik.uni-erlangen.de 



two degrees of freedomji&~— 

In this paper we consider the BC between two gyrating 
charged particles treating the interaction (Coulomb) as a 
perturbation to their helical motions. For electron-heavy 
ion collisions this has been done previously in first-order 
in the ion charge Z for an ion at rest^S and up to Z 2 
for an uniformly moving heavy ion^r— Here we focus 
on second-order perturbation theory and its comparison 
with classical trajectory Monte Carlo (CTMC) simula- 
tions. The present work considerably extends the ear- 
lier studies in Refs. IT3 - il6l where the second-order energy 
transfer was calculated for an electron- heavy ion BC ne- 
glecting the gyration of the ion or assuming identical gy- 
rating particles (e.g. electrons). But in some cases of an 
interaction of magnetized electrons with light ions (e.g. 
protons, antiprotons) at a rather strong magnetic field 
(e.g. in traps) the cyclotron motion of the ion cannot be 
neglected. 

The general expressions for the energy transfers are 
briefly discussed in Sec. HU In Sec. IIIII we discuss the 
perturbative methodology used in this work. Then we 
turn to the explicit calculation of the second-order energy 
transfer for two arbitrary gyrating particle collision and 
without any restriction on the magnetic field. In Sec. IIVI 
the results of the perturbative BC model are compared 
with CTMC simulations. The results are summarized 
and discussed in Sec. |Vj The small velocity limits of the 
energy transfers are derived in the Appendix. 



II. BINARY COLLISION FORMULATION 

For our description of BC we start with the equations 
of motion for two charged particles moving in a homo- 
geneous magnetic field and the related conservation laws 
in general. Next the quantities of interest, the velocity 
transfer, and the energy transfer of particles during the 
binary collision will be introduced and discussed, before 
we turn to the solution of the equations of motion in 
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the subsequent section. For further details we refer to 
Refs. MM- 

We consider two point charges with masses m\ and m 2 
and charges q\e and q 2 e, respectively, moving in a homo- 
geneous magnetic field B = Bh and interacting with the 
potential qiq 2 fi 2 U(r) with ^ 2 = e 2 /47T£ - Here £o is the 
permittivity of the vacuum and r = ri — r 2 is the relative 
coordinate of the colliding particles. 

In plasma applications the bare Coulomb interaction 
t/c( r ) = 1/r is shielded by the surrounding plasma par- 
ticles and the interaction may be modeled by Un(r) = 
e~ r '^/r. Relative velocities which exceed the thermal ve- 
locity lead to an asymmetric interaction potential which 
in general considerably complicates the theoretical de- 
scription. It has been shown, however, that such a dy- 
namic, highly asymmetric interaction potential can be 
replaced with an effective spherically symmetric velocity- 
dependent interaction C/d(i") with a velocity-dependent 
screening length A , 10 i 21 i 22 We adopt these findings for 
our present considerations and assume a spherically sym- 
metric interaction with a given fixed screening length 
in both our analytic expressions and the CTMC simu- 
lations. For the envisaged applications of our given re- 
sults, such as cooling forces, stopping power, etc., which 
typically involves an average over the velocity distribu- 
tion, the screening length has then to be replaced by an 
appropriately chosen velocity-dependent one. 

The quantum uncertainty principle prevents particles 
(for qiq 2 < 0) from falling into the center of these po- 
tentials. In a classical picture this can be achieved by 
regularization of U(r) at the origin. Such a regularized 
potential (pseudopotential) has been derived from quan- 
tum statistical considerations ! 23 ' 24 In our forthcoming 
investigations we take the functional form of this short 
distance correction and, including as well the screening 
contribution, hence use the interaction J7a(r) = (1 — 
e~ T l x )er T l x jr. It should be emphasized, however, that 
the use of this regularized interaction, where A is usually 
related to the (thermal) de Broglie wavelength ; 23 ' 24 here 
essentially represents an alternative implementation of 
the standard (lower) cutoff procedure needed to handle 
the hard collisions in a classical perturbative approach 
where A is taken as a given constant or as a function of 
the classical collision diameter (see Sec. IIV|) . 

Introducing the cyclotron frequencies of the particles 
^cv = \<lv\e-B/m v (with v = 1,2) we start with the clas- 
sical equations of motion 

**(*) ~ Wcv [v„(t) x b] = g ,^-F[r(t)], (1) 

where = \q v \/q V) Qi = 1, g 2 = -1. Here qiq 2 £ 2 F[r(t)} 
(F = —dU/dr) is the force exerted by particle 2 on par- 
ticle 1. In the presence of an external magnetic field, 
the Lagrangian and the corresponding equations of mo- 
tion cannot be separated into parts describing the rel- 
ative motion [r = ri — r2,v = r] and the motion of 
the cm [R = (mirj + m 2 r 2 )/(mi + m 2 ), V = R], in 
gener al & 1 1 ' 1 3 ~ — Introducing the reduced and the total 



masses = 1/mi + l/m 2 , M = mi + m 2 , respec- 
tively, and recalling that v v {t) — V(t) + Qv{n/m v )v{t) 
the equations of the relative and the cm motion are 

v-W3[vxb]=«a[Vxb] + ^-F[r(t)], (2) 

V-wi[Vxb] = jjlj 2 [v x b] . (3) 

The frequencies w\, uj 2 and W3 are expressed in terms 
of the cyclotron frequencies of the particles, io\ — 
(miftWci + m 2 ^ 2 uj c2 )/M, uj 2 = ftw c i - ? 2 w c2 , w 3 = 
(rn 2 <;iU! c i + miq 2 ui c2 )/M . Note that the quantities |wi| 
and \uj 2 \ play the role of the cm and relative cyclotron 
frequencies. The coupled, nonlinear differential equa- 
tions © and d3]) (or Eq. ([1])) completely describe the mo- 
tion of the particles. For solving the scattering problem, 
they have to be integrated numerically for a complete set 
of the initial conditions. 

From Eqs. @ and ^ follow the conservation of the 
parallel component of the cm velocity V(t) ■ b = Voii 
and the total energy but since, in general, the relative 
and center-of-mass motions are coupled the relative E r 
and cm E cm energies are not conserved separately. An 
exception is the case with lo 2 — (or qi/m\ — q 2 /m 2 ) 
where the energies E x and E cm are conserved separately. 

In the general case the rate dE u /dt at which the energy 
E v = m v v 2 /2 of particle v changes during the collision 
with the other particle can be obtained by multiplying 
the equation of motion for particle v by its velocity v„(t). 
The integration of this rate over the whole collision yields 
the energy transfer— _ — 

//•oo 
dkU(k) / [k- V(t)]e lk - rW dt (4) 
J — 00 

assuming that for t — > ±00, r(i) — > 00 and ?7[r(t)] — > 0. 
According to the conservation of total energy we have 
AE 2 — —AEi, as it can be directly seen from Eqs. (fTJ) 
and (HI). In Eq. (Q} the force F(r) has been written using 
a Fourier transformation in space. 



III. PERTURBATIVE TREATMENT 

We now seek an approximate solution of Eqs. (P)-© 
by assuming the interaction force between the particles as 
a perturbation to their free helical motion. For the case 
of electron-ion (without cyclotron motion of the ion) and 
electron-electron scattering the corresponding consider- 
ations and derivations are discussed in Refs. [IH and [iH 
We therefore focus here on the general case. 

We look for a solution of Eq. (JTJ) for the variables 
iV,v„ = r v , or alternatively, of Eqs. © and for 
the variables r, v = r and R, V = R, in a perturbative 
manneri^, i.e. r„ = r^ 0) +r^ 1) + ... andR = R +Ri+..., 
r = ro + r^ 1 ' + ... Starting point is the zero-order unper- 
turbed helical motion of two particles in the laboratory 
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frame with vi. -* (t) = iu'' if) and 

r( / °)(t) = R„ + b^||t (5) 
+a u [u„ sin (u cu t) + ^[bx u„] cos {ui cv t)\ , 

where u„ = (costly, sin tp„) {ip v is the initial phase of 
the particle v) is the unit vector perpendicular to the 
magnetic field, v§ v u and v§ v \_Vl v (with v§ v \_ 0) are 
the unperturbed velocity components parallel and per- 
pendicular to b, respectively. Here a„ = vq v ±/oj cll is the 
cyclotron radius of the particle v. It should be noted that 
in Eq. ((5]) the variables and R„ are independent and 
are defined by the initial conditions. The unperturbed 
cm (Ro(i),Vo(£) = Ro(i)) an d relative (r (i),v (i) — 
i"o(t)) coordinates and velocities can be easily found from 
Eq. ([5]). In general, the cm and relative coordinates and 
velocities involve two harmonic oscillations with differ- 
ent frequencies, amplitudes and phases. Therefore in the 
plane perpendicular to b these quantities cannot be rep- 
resented in the form of a simple cyclotron motion with 
constant amplitudes and phases as in the case of electron- 
electron (or electron-positron) collisions^ 

The equation for the first-order velocity correction is 
given by 



(0), 



vW(i)- W ,[vW(t)xb] = fo 



-F[r (t)], (6) 



into Eq. ([7]). This quantity is then averaged with re- 
spect to the initial phases of the particles tpi and ip 2 
and the azimuthal angle <& s of the impact parameter s. 
The obtained angular integrals are evaluated using the 
Fourier series of the exponential function. After averag- 

(2) 

ing the energy transfer AE{ with respect to ipi and ip 2 
the remaining part will depend on <5(fcj| + fcj|), i.e. the 
component of k + k' along the magnetic field b. This 5- 
function enforces k + k' to lie in the plane transverse to 
b so that e i ( k+k ')- R -(5(fc|| +&[,) = e iQ ' s <5(A;|| + fc[,), where 
Q = kj^ + k^ . The result of the angular averaging finally 
reads 

{AE (2) ) = Ttiq{q& I dkdk/u (k) u (k , } Jo {Qs) 



H \v r \\\ 



l -S')( r "Ji+ m ft) 



n,m— — oo 



x ( fenVoii + —nuj. 
\ m 2 



fx \ I 2fc 2 i?„. m (fc^, fcj_) 



ci muj c2 



k' x k A 



it) 



\Cn,m ^0) 

2iAH njrn (k x ,k' x )sm(0-e') 



(i 



(k ± ,k' ± ) 



m 2 uj c2 



(8) 



where the solutions can be given by an integral involv- 



ing the force F[r (t)] with r (i) = (t) - r 2 v '(t) and 
the unperturbed trajectory ([5]) by similar expressions as 
Eqs. (43)-(46) of Ref. (see also Ref. GJ and[H). 

The first- (A£^ 1} ) and the second-order {Ae[ 2) ) en- 
ergy transfers of the particle I can be evaluated using 
general Eq. Q . Here we consider only the energy change 

AE^ since the angular averaged Ae[^ vanishes due to 
the symmetry reasons*^— We obtain 

AE {2) = qi q 2 f f dkU(k) f e lkro(t) di 



x{(k.rW(t)][k.Vo(i)]-i[k.V 1 (*)]}. 



(7) 



Let us recall that here r^'(t) and Vi(t) are the first- 
order relative coordinate and the cm velocity corrections, 
respectively. 

We now introduce Jhe variable s = R r ^ which is the 
component of R r = Ri — R2 perpendicular to the rela- 
tive velocity of the guiding centers of two particles u r yb, 
where v r \\ = i>oi|| — ^02||- From Eq. ([5]) we can see that 
s is the distance of closest approach for the guiding cen- 
ters of the two particles' helical motion. For practical 
applications the energy change is given by the average 
of AEi with respect to the initial phases of the particles 
if 1 and ip 2 and the azimuthal angle i9 s of s. This aver- 
aged quantity is denoted by (AEi) in the forthcoming 
considerations. 

(2) 

To evaluate the second-order energy transfer AE{ 
we have to insert Eq. §5$ and the solution of Eq. © 



Here H n<m = J n (kj_cii) J n (k' ± cii) J m (kj_a 2 ) J rn (k' ± a 2 ) and 

^n,m — ^*n— l,m -^n+l.m; T^n^m — -^n,m — 1 -^n,m+l; 

A = fJ>Si/miu c i + liQz/rnzw&t where J n are the Bcsscl 
functions of the nth order. Also Cn,m — nuj ci + muj c2 + 
k\\v r \\ , fc|| = k b and k± are the components of k parallel 
and transverse to b, respectively, tan (9 = k y /k x . The se- 
ries representation (|8]) of the second-order energy transfer 
is valid for any strength of the magnetic field. 

For most applications it is also useful to integrate the 
averaged energy transfer, (AE\ ), with respect to the 
impact parameters s in the full two-dimensional (2D) 
space. We thus introduce an energy loss cross section 
(ELCS)&±&±i-.i£ through the relation 



(AE[ 2) )sds. 



(9) 



As a results from the s-integration of the energy trans- 
fer ([5]) one obtains an expression for a which represents 
an infinite sum over Bessel functions. Moreover, assum- 
ing regularized interaction and performing k integration 
in a yields an infinite sum over modified Bessel functions 
(see, e.g., an example for ion-electron collision in Ref. I16T) . 
For arbitrary axially symmetric interaction potential sim- 
ilar expression is derived in the Appendix [see Eq. (IAI[) ]. 
However, for practical applications it is much more con- 
venient to use an equivalent integral representation of the 
ELCS which does not involve any special function. This 
expression can be derived from the Bcsscl-function rep- 
resentation of a using the integral representation of the 
Dirac 5 function as well as the summation formula for 
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T n e mip Jn( a )^ The energy transfer a after lengthy but 
straightforward calculations then reads 

a = [ ^a(tp), a(tp) = a\\{tp) + a±(tp), (10) 
Jo Z7r 



with 



CT|| (<p) = - 



(2tt) 2 g gg|^ 4 A 2 Vb|| 

^11 



f tdt J dk|C/(k)| 2 



x [k\ + k\(j> (t)]k\\ sin (fc||At) J [2fc^ft (t)] , (11) 



OU. (V) 



miTO2'i 



2 ; 2 

X 



Ji [2fc ± ft (£)] 



x cos (fc|, A<) <^ A 2 Gi (t) [fcf + fci0 (t)] 



2A;_l71 (i) 



+G 2 (t) J [2M(t)] 



(12) 



where 5„ = |t> r ii|/a; el , are the relative pitches of the par- 
ticles helices, divided by 2ir and 77„ = \/8 v = Lo cv X/\v r \\\. 
In Eqs. (fTTj) and (fl~2j) the time i is scaled in units A/|t> r |||, 
where the length A is specified in the next section. In 
addition in Eq. (1101) the energy transfer a(tp) has been 
split into two parts which correspond to the cm motion 
along [<7|| (ip)] and transverse [^_l(<^)] to the magnetic 
field, where tp = tp± — tp 2 is the difference of the initial 
phases of the particles. Also 



0(t) = 
TZ 2 (t) 



/i sin (771 1) fx sin (772*) 



mi 



2 »7l* 
:n 

2 



m 2 



- a Q sin 



2 ^ 
2 



. mt ■ mt 

-la\a 2 sin — - sm cos 9?, 



(13) 



(14) 



^ /.X 2 Sm (Vlt) 2 

GiW=™i%i t 

771 rj 



sin (772 t) 
TO 2f 2J_ ; H 2v i ± v 02 ± 



f/2* 



sin^ 771 i 

x cos tp m 2 cos mi — cos — 

\ 771* 2 772* 2 



sm 2^ 772 < 



G 2 (t) = /ju r || 



2 cos (772*) — cos (771 1) 



t 2 



(15) 
(16) 



Here the quantity lZ(t) represents an effective relative cy- 
clotron radius of the two particles. The second term in 
Eq. (TT2"j) proportional to the function G 2 (i) is the con- 
tribution of the cm velocity perturbation, see the second 
term of Eq. ([7]). This perturbation and hence this corre- 
sponding term was absent in our previous considerations 
for electron-electron and electron-ion collisions. 

An expression similar to Eq. ([5J (or Eqs. (fLOjl - (fT2|l ) 
has been obtained for electron-electron collision 15 and for 
electron-heavy ion (without cyclotron motion of the ion) 
collision*^ - — In contrast to these specific cases Eqs. @ 
and (fK>)) -(fl"2" j) are valid for the collision of two arbitrary 
charged particles. The results obtained in Refs. [l4T - [l6l 



can be easily derived from these general expressions. The 
electron-electron case is recovered assuming that m\ = 
m 2 = m, qx = q 2 = q, i.e., lo c1 = u> c2 = lo c , T)i = r/ 2 = 77, 
fi = m/2. Then the term proportional to the function 
G2 (t) in Eq. (1121) vanishes and the remaining expressions 
coincide with the results of Ref. [H. It should be noted 
that for identical particles (e.g. electrons) the effective 
cyclotron radius is given by lZ(t) = a sia(T]t/2), where 
a 2 = a\ + a 2 — 2a\a 2 cos tp is the cyclotron radius of the 
particles in the relative framed For electron-heavy ion 
collision we assume that the ion mass tends to infinity, 
?7i2 — > 00, and vq 2 ± = (ion moves along the magnetic 
field direction). Therefore /1 — > toi = to, lu c2 — > (772 — >• 
0). In this case the transverse ELCS vanishes and the 
longitudinal cross section coincides with the results of 
Ref. OJ. 

Next we also consider the ELCS W\\(tp) and W±(tp) for 
vanishing cyclotron radii, a\,a 2 — > 0, i.e. when ini- 
tially the particles move along the magnetic field (vqi± = 
tjo2± = 0). In this limit <tj_ (tp) 7^ as in the cases 
of two identical particles collision or ion-electron colli- 
sion. For any axially symmetric interaction potential, 
t/(k) = U(\kn I, kx), the straightforward integration in 
Eqs. (HQ) and flT2]) yields 



CT || O) = - 



2g^%i 

V r\\ 



TOl 



TO 2 



(17) 



^) = ^-[M^)-M^)]> ( 18 ) 



where 



Jb (k) 



(27T) 4 



U 2 (K,k ± )k 3 j_dk ± . 



(19) 



The averaged energy transfer, Eq. ©, can be evaluated 
without further approximation for any axially symmetric 
interaction potential. In this case the energy transfer can 
be represented as the sum of all cyclotron harmonics as 
it has been done for ion-electron interaction in Ref. [l6| 
and for electron-electron interaction in Ref. IT5I 

In the following we consider the regularized screened 
potential introduced in Sec. ITTI with 



C/ R (fc||,fcj 



1 



1 



(2tt) 2 \kl + n 2 kl+ X 2 J ' 



(20) 



where k 2 = kj + A -2 , X 2 = *f + d' 1 = A" 1 + A . 

As we discussed above the energy transfer (JSJ) must be 
integrated with respect to the impact parameters s for 
practical applications. For general interaction potential 
this is given by Eqs. ((9|)- (fT2|) . In general for a study of 
the convergence of the s-integrated energy transfers we 
note that the case with some value s = s c is most critical 
for the convergence of the ELCS. This is intuitively clear 
as the gyrating particles at \a± — a 2 \ < s < a± + a 2 may 
hit each other on such a trajectory (see Ref. [l6| for some 
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explicit examples). This should not matter for the po- 
tential (|20p . which has been regularized near the origin 
for exactly that purpose. 

For the present case of the regularized interaction po- 
tential, substituting Eq. (g(| into Eqs. (JTU) and ([12]). we 
obtain 



+ 



x 2 - 1 
4x 2 

x 2 -l 



T 2 [R(t),t,t] 



Ji[xii(t),jrt,t] 



J" 2 [x#(t),xi,t] 



qhU 4 



(21) 



4 



G x (t) 



x _F 3 [x.R(*),x£,t] 
+G 2 (t)/e- B W 



x 2 - 1 



B?{t)\ 
T 4 [R(t),t,t 



,-R(t) 



4x 2 

X 2 - 1 

J" 5 [i?(t),t] H 



•F 4 [x.R(i),xt,t] 
4 



x 2 - 1 



T 6 [R(t),t] 



(22) 



<il(t) 



;F 5 [xE(f),xi] 



4x 2 
x 2 - 1 



J" 6 [x.R(t),xi] 



Here i? 2 (i) = t 2 + (4/A 2 )ft 2 (t), x = X/d = 1 + A/A, and 



Ji(J2, C, *) - 2 + 2i? - R 2 + [1 - 0(t)] 



C 2 

^r(R 2 + 3R + 3) 



R 2 + R + 1 



i? 2 



j- 2 (i?,C,t) = i?.+ i-^[i-0(t)] 

-^(R 3 + 6R 2 + 15R+ 15) 
fir 



J- 3 (i?, C, *) - 2 + 2i? - i? 2 + [1 - 0(i)] 



(23) 



J? + 4R 2 + 9R + 9 



(24) 



R 2 -R-l 



i? 2 



(i? 2 + 3i? + 3) 



74(i2,c,t) = -R + i--^[i-#*)] 



c 2 



(R 3 + 6R 2 + 15R+ 15) 



I? 2 

J- 5 (i?,C) =i? 2 



l-i?+-^(i?+l) 



^(i?, C) = i? 2 + i? + 1 - 7^(i? 2 + 3i? + 3). 



i? 2 ' 



(25) 

R 3 + 2R 2 + 3R + 3 
(26) 
(27) 
(28) 



Next we consider the ELCS cr\\(<p) and (r±(ip) for van- 
ishing cyclotron radius, a\, a 2 — >■ 0, i.e. when initially the 
particles move along the magnetic field (i>oi_l = ^02± = 



0) which are given by Eqs. (|17l) - (fH))) . Substitution of 
Eq. (J20I) into Eq. (Ql) yields 



T (k) 



2k 2 A 2 + x 2 + 1 k 2 A 2 + x 2 

■ In — ;ttt; 1 • 



2(x 2 - 1) 



K 2 A 2 + 1 



(29) 



Thus the ELCS at at, a 2 are then given by Eqs. (|T7|) 
and (HHJ), where the function J-o(k) is defined by Eq. (|29D . 

Equations (fTT)) and ([TH| with Eq. ((29)) are approxi- 
mately valid also for finite cyclotron radii a\ and a 2 , as- 
suming that the longitudinal velocity i? r ii is larger than 
the transversal ones, i>oi± and vq 2 ±. Indeed, in this case 
R{t) ~ i since 5 X , <5 2 > ai,a 2 and in Eqs. (JST) and ((22)) 
the transversal ELCS ct_l(<^) can be neglected compared 
to the longitudinal one. This indicates that in the high 
velocity limit with v r \\ 3> Woi_L) ^02± the transversal mo- 
tion of the particles as well as its cm transversal motion 
are not important and can be neglected. Since only the 
contribution of small t is important the function R(t) is 
approximated by R(t) ~ t. Using this result for R(t) 
from Eq. $Z2§ and Eqs. $T7$ and (QJ§ with $Z§§ in the 
high velocity limit we obtain within the leading term ap- 
proximation for the ELCS 

_ ( v 2q 2 q 2 ^V 0ll 

a \\( ( P)- — A (x), (30) 



[QiA(x) + Q a Ai(x)], (31) 



where Q 2 = (A 2 /2)(w 2 x - w 2 2 ), Qi = - 



(m 2 /n)v 2 2± 



((m 2 - mi)/^i)voi±«o2± cosy, 
1 



A(u 



1 



In u — 1 , 



1 



Ai (u) = 1 + -=■ - 



■ In 'i 



(32) 
(33) 



Note that W\\(<p) is isotropic, i.e., do not depend on (p 
while a±_{ip) contains a term which is proportional to 
cosy. Also, in the high velocity limit the parallel ELCS, 
ctm (</?), does not depend on the magnetic field strength. 



The ELCS decay as a»(ip) 



and cr±(ip) 



Finally we briefly turn to the case of small relative ve- 
locity, w r ii <C Uoi_l, i>02X- It should be emphasized that 
Eqs. (12"T)) and f2"2")) . are not adopted for evaluation of the 
ELCS at small velocities. For this purpose it is conve- 
nient to use an alternative Bessel-function representation 
of the ELCS as shown in the Appendix. In addition, it 
is expected that the limit of small u r ii is the most criti- 
cal regime for a violation of the perturbation theory em- 
ployed here. Therefore explicit analytical expressions in 
this limit can be useful for an improvement of the per- 
turbation theory by comparing the analytical results with 
numerical simulations, see Sec. IIVI 



IV. RESULTS. COMPARISON WITH SIMULATIONS 

A fully numerical treatment is required for applica- 
tions beyond the perturbative regime and for checking 
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FIG. 1. (Color online) Top panels, the ELCS cry for electron-ion (Z = 6 and m^/mx ~ 2.2 x 10 4 ) collision at Vi±/v s = 10~ 3 
(left panel) and Vi±/v s = 3.0 (right panel). The curves with and without symbols correspond to CTMC simulations and the 
second-order perturbative treatment, respectively. Also A/a se = 45, v e ±/v s = 0.5625 (solid lines), v e ±/v s = 2.812 (dashed 
lines), v e ±/v s — 7.987 (dotted lines). Bottom panels, same as in top but for a±. 



the validity of the perturbative approach outlined above. 
In the present case of binary collisions of two arbitrary 
particles in a magnetic field and with the effective in- 
teraction Un(r) ([20)1 the numerical evaluation of the BC 
energy loss is very complicated, but can be successfully 
investigated by classical trajectory Monte Carlo (CTMC) 
simulations^— In the CTMC method^ the trajectories 
for the relative motion between the particles are calcu- 
lated by a numerical integration of the equations of mo- 
tion (0) and ([3]), starting with initial conditions for the 
parallel v r \\ and the transverse v ± = v ij_ — v 2± rel- 
ative velocity. The required accuracy is achieved by us- 
ing a modified velocity- Verlet algorithm which has been 
specifically designed for particle propagation in a (strong) 
magnetic field, 27 >M and by adapting continuously the 
actual time-step by monitoring the constant of motion 
E = E\+ E<2, = E T + E cm . The resulting relative devia- 
tions of E are of the order of 10~ 6 — 10~ 5 . 

The desired average over the initial phases and the 
impact parameter s is performed by a Monte Carlo 



samplin g 29 i 30 of a large number of trajectories with dif- 
ferent initial values. The actual number of computed tra- 
jectories (typically 10 5 — 10 6 trajectories) is adjusted by 
monitoring the convergence of the averagi ng p rocedure. 
For further details we refer to Refs. [TH andllq. 

For the forthcoming discussion we put the equation 
of the relative motion of two particles in a more appro- 
priate dimensionless form by scaling lengths in units of 
the screening length A and velocities in units of a char- 
acteristic velocity v s defined by v 2 — \q1q2\f 2 / This 
velocity gives a measure for the strength of the Coulomb 
interaction with respect to the (initial) kinetic energy of 
relative motion \jLv 2 j1. For v r < v s the kinetic energy 
is small compared to the characteristic potential energy 
<M2|^ 2 /A in a screened Coulomb potential and we ex- 
pect to be in a nonperturbative regime. A perturbative 
treatment on the other hand should be applicable for 
v r » v s . 

The scaled version of Eqs. (P)-© depends on the four 
dimensionless parameters a s \/X, a S 2/A, 7712/7711 and A/A, 
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and the initial conditions with positions scaled in A and 
velocities in v s . Here a s \ = v s /ui c i and a S 2 — v s /uj C 2 are 
the cyclotron radii for v±± = t>2± = v s and the parame- 
ters a s i/X oc a S 2/X oc v s /B represents a measure for the 
strength of the magnetic field compared to the strength 
of the Coulomb interaction [which is oc v 2 s ]. The ratio 
A/ A describes the amount of softening of the screened 
interaction at r — > with qxq^^U^ir — > 0) — > q^fc 2 /X. 

In the analytical perturbative approach we thus apply 
the same scaling of length and velocities and introduce 
for two particles collisions the dimensionless ELCS 

C-a = \2 [ 7T Wa (v) . (3 4 ) 

where a =\\, _L, Vii = Vqii and V± =v s . 

Next we specify the cutoff parameter A which is a mea- 
sure of softening of the interaction potential at short dis- 
tances. As we discussed in previous section the regular- 
ization in the potential (|20p is sufficient to guarantee the 
existence of the s-intcgrated energy transfers, but there 
remains the problem of treating hard collisions. For a 
perturbation treatment the change in relative velocity 



must be small compared to v r and this condition is in- 
creasingly difficult to fulfill in the regime v r —> 0. This 
suggests a physically reasonable procedure: the poten- 
tial must be softened near the origin. In fact the pa- 
rameter A should be related to the de Broglie wavelength 
which is inversely proportional to v r . Here within classi- 
cal picture we employ in a perturbative treatment the dy- 
namical cutoff parameter x(w r ii) = 1 + X/ X(v r \\)^ where 
A 2 (w,.||) = Cbl{v A ) + A^ with &o(tV||) = Av s 7(^n + «o), 
v o = u oij_ + V Q2±- H ere Ao is some constant cutoff pa- 
rameter, and &o( w r|| ) is the distance of closest approach of 
two charged particles in the absence of a magnetic field. 
Also in A(tV||) we have introduced a fitting parameter 
C ~ 0.292. In Ref. 15 this parameter is deduced from 
the comparison of the ELCS © with an exact asymptotic 
expression derived in Ref. |3l| for the Yukawa- type (i.e., 
with A — > 0) interaction potential. As has been shown 
in Ref. [TH the second-order ELCS for electron-electron 
and electron-ion (but wihtout gyration of the ion) colli- 
sions with dynamical cutoff parameter A(u r ii) excellently 
agrees with CTMC simulations at high velocities. The 
CTMC simulations have been carried out with constant 
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FIG. 3. (Color online) Top panels, the ELCS cxii for electron-ion [Z — 6 and mj/mi ~ 2.2 x 10 4 ) collision at Vi±/v s = 0.5 (left 
panel) and Vi±/v 3 — 3.0 (right panel). The curves with and without symbols correspond to the CTMC simulations with and 
without ion gyration effects, respectively. Also X/a se — 225, v e ±/v s — 0.5625 (solid lines), v e ±/v s — 2.812 (dashed lines) and 
v e ±/v a — 7.988 (dotted lines). Bottom panels, same as in top but for a±. 



X = %o <C A, that is, the interaction is almost Coulomb 
at short distances. 

The ELCS are presented in Figs. HE These results are 
obtained for the BC of an electron (particle 1, q\ = — 1) 
with an ion (particle 2, q<i = Z) in a strong magnetic 
field. Shown are <tii and a± as functions of v r u / v s for fixed 
transverse velocity of the ion vq2± = Vi± and the strength 
of the magnetic field \/a se = B/B s (with a se = a s \, 
B s = miV s /eX) and varying the transversal velocity of 
the electron vqi± = v e ±. For each triplet of fixed Vi±, 
v e j_ and B the cyclotron radii ai, ai and a S 2 = a S i can 
be easily determined. Both the CTMC and second-order 
calculations have been done for a regularized potential 
[/r with vo = Xq/X = 10~ 4 . Note that in all cases shown 
in Figs. HE] we have a se /X <C 1 and the ELCS is not 
sensitive to this parameter when a se /X — > (see, e.g., 
Fig.®. 

Comparisons of the ELCS determined by the CTMC 
simulations and the second-order perturbative treatment 
Eqs. (|2l ]) . (|22 )l . and (|33 | are presented in Figs. Q] and H 
It is clearly observed that in the regimes of large rela- 



tive velocities the second order perturbative treatment 
agrees almost perfectly (i.e. within the unavoidable nu- 
merical fluctuations) with the CTMC results. In addi- 
tion in the limit of very large velocities v r u/v s 3> 1 the 
ELCS (Tii calculated either within perturbation theory 
or CTMC method with different strength of the mag- 
netic field and Vi±, v e ± converge to the same value. This 
behavior agrees with the predictions of the asymptotic 
Eq. pop which is independent on B and Vi±, v e ±. Also 
the smaller the transversal velocities, the better is the 
convergence to the regime of Eq. (|30[) . At large rela- 
tive velocities the CTMC and second-order <r± shown 
in Figs. Q] and [5] agree with Eq. with the asymp- 
totic Eq. (|3ip . Since at v r \\/v s 3> 1 the quantity a± 
behaves as a± ~ vfj_ for fixed magnetic field and v e ± it 
will increase with the transverse velocity as shown 
in Figs. [1] and [H At small velocities with v r \\/v s < 8 
the second-order treatment considerably deviates from 
CTMC simulations. Here the second-order ELCS are 
given by approximate expressions (|A6p and (|A7p where 
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FIG. 4. (Color online) The ELCS cm for electron-antiproton {Z = — 1 and mj/mi = 1836) collision at Vi±/v s — 2.0 (left 
panel) and Vi±/v s — 3.0 (right panel). The curves with symbols correspond to CTMC simulations. Also X/a se = 2.2 x 10 4 , 
v e ±/v s — 0.176 (filled circles), v e ±/v a — 2.499 (squares), v e ±/v s = 8.8 (open circles). The dashed curves are obtained employing 
a model suggested in Ref. I la . 



the parameter yt at small relative velocities is given by 
>t ~ x(0) = 1 + A/A(0) and A(0) is the dynamical cutoff 
A(w r ||) at v r ii = 0. Note that at finite cyclotron radii of 
the particles the quantity A(0) is a constant depending 
on the value of vq and the transversal velocities. How- 
ever, for vanishing cyclotron radii (as, e.g. in Eqs. (|17[) 
and (fTS| ) the cutoff parameter at small velocities behaves 
as A/A ~ (v s /v r \\) 2 . The quantity (x — l) 2 involved in 
Eqs. flMJ and (|A"7j) falls as ~ {v r \\/v s ) A . This results in 
a strong self-cutting at small velocities. Thus employing 
the cutoff A(iV||) the second order ELCS an is strongly 
reduced and decreases as <j\\ ~ tr.i and cr|| ~ ln(l/w r ||) 
at ai = <i2 = and a\,a2 ^ 0, respectively. The second 
term of a± in Eq. (IA7|) does not contain a term (x — l) 2 
and diverges as a j_ ~ v~ 2 , see Figs. [1] and [5] In this small 
velocity regime the second-order perturbative treatment 
is clearly invalid and a nonperturbative description is re- 
quired. 

To highlight the importance of the ion gyration in the 
presence of very strong magnetic field we demonstrate in 
Fig. H the ELCS obtained with CTMC simulations with 
and without ion gyration. That is, in the latter case we 
assume that ui2 — > oo and the ion moves with rectilinear 
trajectory with the same velocity component Vi± as in 
the case of ion gyration. It is seen that the ion gyration 
is important at small v r \\ and the discrepancy between 
two approaches increases with Vi±. 

In Fig. 0] we compared the CTMC results for a\\ with 
the model (dashed lines) given in Ref. [TH for a repulsive 
(Z < 0) ion-electron interactions. Here the ions (antipro- 
tons) are strongly magnetized with a^/A ~ 0.17 (left 
panel) and a-i /A ~ 0.25 (right panel). Assuming that 
the magnetic field is infinitely strong this model com- 
pletely ignores the cyclotron motion of the particles and 
they move along b. It has been showr*i£ that for re- 



pulsive interaction the magnetic field together with the 
interaction potential forms a potential barrier because of 
the particles motion is effectively one-dimensional. In 
this case the relative velocity transfer is Awii = — 2u r ii 
which corresponds to a reversion of the initial motion, 
i.e. to a backscattering event. Then the energy transfer 
is AE 1 = -2/iVo|| v r \\Q(vl - v^)Q(s m - s), where Q(z) 

is the Heavyside function, v 2 c = 2\q\q2\^ 2 / ^iX and s m is 
determined from equation \qiq2\^ 2 Un(s m ) = fiv 2 ^/2. It 
is seen that in Fig. @] the agreement with CTMC simu- 
lations is quite satisfactory even for finite cyclotron radii 
and magnetic field. However, with increasing Vi± the 
CTMC simulations show a more involved picture as, e.g., 
in Fig. right panel, than the predictions of this sim- 
ple model. In the CTMC simulations the ELCS shrinks 
strongly at v r u > v s with increasing relative velocity and 
Vi± , i"e_L • At strong but finite magnetic field the hard col- 
lisions like backscattering events may also occur but the 
transverse dynamics of the particles will reduce the do- 
main of the backscattering events^. With increasing v r \\ 
this domain will be further shrunk and finally the scatter- 
ing may occur only in the regime where a strong magnetic 
field may strongly reduce the energy transfer J^—. 



V. CONCLUSION 

In this paper we have investigated the binary collisions 
(BC) of two gyrating charged particles in the presence 
of constant magnetic field employing second-order per- 
turbation theory and classical trajectory Monte Carlo 
(CTMC) simulations. A case with strongly asymmet- 
ric masses of the particles have been considered in de- 
tail. The second-order energy transfers for two-particles 
collision is calculated with the help of an improved BC 
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treatment which is valid for any strength of the magnetic 
field and involves all cyclotron harmonics of the particles 
motion. For further applications (e.g., in cooling of ion 
beams, transport phenomena in magnetized plasmas) the 
actual calculations of the energy transfers have been done 
with a screened interaction potential which is regularized 
at the origin. The use of that potential can be viewed as 
an alternative to the standard cutoff procedure. 

For checking the validity of the perturbative approach 
and also for applications beyond the perturbative regime 
we have employed numerical CTMC simulations. These 
CTMC calculations have been performed for a very 
strong magnetic field and in a wide range of v r u and for 
a small regularization parameter, that is, for an interac- 
tion which is rather close to Coulomb at short distances. 
Within the second-order treatment we have introduced a 
dynamic cutoff parameter which substantially improves 
the agreement of the theory with CTMC simulations. 
From a comparison with the nonperturbative CTMC sim- 
ulations we have found as a quite general rule which is 
widely independent of the magnetic field strength that 
the predictions of the second-order perturbative treat- 
ment are very accurate for v r ii/v s > 8 for all studied pa- 
rameters and cases. In contrast, for low relative velocities 
v r\\/ v s ^$ 8 the results obtained from perturbation theory 
strongly deviate from the CTMC simulations. We have 
also tested the exact analytical model derived for repul- 
sive interaction and an infinitely strong magnetic field in 
Ref . flU bv comparing it in Fig. [J] with the CTMC simula- 
tions and found that the agreement is rather satisfactory 
even for finite (but strong) magnetic fields. 

We believe that our theoretical findings will be use- 
ful for the interpretation of experimental investigations. 
Here, it is of particular interest to study some macro- 
scopic physical quantities on the basis of the presented 
theoretical model such as cooling forces in storage rings 
and traps, stopping power of ion beams as well as trans- 
port coefficients in strongly magnetized plasmas. These 
studies require an average of the energy or velocity trans- 
fers with respect to the velocity distribution of the elec- 
trons. The cooling forces obtained by the perturbative 
approach are expected to be quite accurate if the low 
velocity regime only slightly contributes to the v r aver- 
age over (AE). That is, if the typical v r n, given by the 
maximum of the thermal electron velocity and the ion 
velocity, are large compared to v s , as it is usually the 
case for, e.g., electron cooling in storage rings. 



Appendix A: The energy transfer in a small velocity limit 

For the second-order BC treatment the most critical 
situation is the small velocity regime where we expect 
some deviations from the nonperturbative CTMC simu- 
lations. For the improvement of the theoretical approach 
it is therefore imperative to investigate the energy trans- 
fer in the small velocity limit, \v r \\\ <C foi_Li ^02^, or al- 
ternatively Si, 82 <C a\, 0,2- In principle this limit can be 
evaluated using the integral representation of the ELCS, 
Eqs. (|21[) and (|22j) . However, while these expressions are 
very convenient to calculate the high velocity limit of the 
energy transfers (see Sec. IIII[) they are not adopted for 
the evaluation of the small velocity limit due to the os- 
cillatory nature of the function R(t) at v r \\ — s- 0. In this 
Appendix we consider instead an alternative but equiva- 
lent expression for the ELCS. For the axially symmetric 
interaction potential the ELCS can be evaluated using 
Eqs. © and ©. We refer the reader to Refs. for 
details. The integration of Eq. ([8]) with respect to the im- 
pact parameter s yields the two-dimensional 5 function, 
<5(kj_ +kj_), which combining with 5(k\\ + fcy) in Eq. ([5]) 
yields a three-dimensional S function 8(k + k'). The k' 
integration in the energy transfer can be then performed 
exactly. Furthermore it can be shown that the ELCS 
is determined by the imaginary part of the integrand in 
Eq. (J5J), which is expressed by the functions (5[Cn, m (k)] 
[see, e.g., Eq. (57) of Ref. [3 for ion-electron collision]. 
This allows to perform the fen integration. The final re- 
sult reads 
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where n nm = £l nm /\v r \\\ = n/Si + m/5 2 , Q n m = nu) c i + 
mw C 2, and the function G> nm (k\\) is defined as 

$nm (fc||) = ^~r- / C/ 2 (fc||,fc ± )j2(fc_ L ai)J ? 2 „(fc i a 2 )fc_ L dfcj 
4 Jo 

(A3) 

Equations (jAll) and (IA2|) are equivalent to the integral 
representations (fTUjl - ffT^ll of the ELCS, respectively. Also 
in Eqs. (|A1|) and (IA2p the function ( & rim (fc||) is taken at 

&II — ^nm- 
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For evaluation of the ELCS mi and a± at small relative 
velocities v r u we note that the terms with n nrn — [in 
particular the term with n = m = 0] do not contribute 
to the parallel ELCS, cry, and they must be excluded 
from the summation in Eq. (|A1|) . Similarly the term 
with n = m = must be excluded from the summation 
in ELCS <7j_. Equation (| A2|) can be split into two parts 
with K nm ^ and n nm — 0. Therefore in Eq. (| Al|) and in 
the first part of Eq. (|A2[> (with K nm ^ 0) in the limit of 
small velocities the quantity n nm becomes very large. For 
the regularized interaction potential (|20|) from Eq. (|A3|) 
at fcyai, fen a.2 S> 1 in the leading order we obtain 



(^-l) 2 



17 

ln (l fc |||A) ~ J^+Snm 



(A4) 



Here u± = a±/X, u 2 = a 2 /A, 

7r 2 ui u 2 (ui x) (u 2 x) 



(x 2 + l) 2 J XdX - 
(A5) 

Substituting the expression (lA4|) into Eqs. (fA"Tj) and (|A"2|l 
in the lowest order with respect to u r ii we arrive at 
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In Eqs. (IA6|) and (|A7[) the terms with fi„ m = in the 
n. m summations must be excluded. The second term in 
Eq. (1A7[) proportional to is the contribution of the 
terms with n nm = 0. Note that this relation requires that 
w C 2/w c i = p/q is the rational fraction (p, q are two arbi- 
trary integers) and to = —kq, n — kp with k = ±1, ±2, ... 
If w c2 /u; c i is not rational fraction the second term in 
Eq. (|A7I) must be omitted. Thus, the small-velocity 
limit of the ELCS strongly depends on the parameter 
w c2 /w c i = Iffe/siK^i/m^). In addition, Eqs. (|A6I) and 
(|A7[) are not valid at w c i = ui c2 , i.e. in the case of two 
symmetrical (e.g., electron-electron) or antisymmetrical 
(e.g., electron-positron) particles collisions. In particular, 
comparing Eqs. (|A6|) and (jA7[) with similar expressions 
derived for electron-electron collisions^ we then obtain 
that at small relative velocities the ELCS behave here 
as cry ~ w 3 n ln(l/|u r |||) and a ± ~ ln(l/|u r |||) + C/v 2 r]] 
(where C is some constant) while in the electron-electron 



case cr ee || ~ w 2 || and <7 ee ± ~ + ^W^y. This spe " 
cial case with w c i = w c2 = w c can be evaluated start- 
ing from Eqs. (|Aip and (|A2|) . Since n nm — (n + m)/5 
(where S — |u r ii|/a; c ) we introduce a new summation 
variable n — » n — to. Then the to summation can 
be performed exactly using the summation formula for 
Y,m mlj n- m (a)Ji(b) with £ = 0,1, 2. 25 Similarly one 
can evaluate the second term of Eq. (|A7|) with p = q = 1 . 

Finally, we consider the small velocity limits of au (tp) 
and ax(ip) at vanishing cyclotron radii, a\ = a 2 = [see 
Eqs. (Q3| and (QSJ) with Eq. (|29j) ]. 
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In this case cry ((f) 



and ctj^ (</?) 



at small rela- 



tive velocity v r || [cf. with Eqs. (|A6I) and (|A7[)]. 
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